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Abstract The behavior of neurons can be modeled by the FitzHugh-Nagumo oscillator model, consisting of 
two nonlinear differential equations, which simulates the behavior of nerve impulse conduction through the 
neuronal membrane. In this work, we numerically study the dynamical behavior of two coupled FitzHugh- 
Nagumo oscillators. We consider unidirectional and bidirectional couplings, for which Lyapunov and isope- 
riodic diagrams were constructed calculating the Lyapunov exponents and the number of the local maxima 
of a variable in one period interval of the time-series, respectively. By numerical continuation method the 
bifurcation curves are also obtained for both couplings. The dynamics of the networks here investigated are 
presented in terms of the variation between the coupling strength of the oscillators and other parameters 
of the system. For the network of two oscillators unidirectionally coupled, the results show the existence 
of Arnold tongues, self-organized sequentially in a branch of a Stern-Brocot tree and by the bifurcation 
curves it became evident the connection between these Arnold tongues with other periodic structures in 
Lyapunov diagrams. That system also present multistability shown in the planes of the basin of attractions. 

Key words. FitzHugh-Nagumo oscillator - Periodic structures - Lyapunov diagrams - Isoperiodic dia¬ 
grams - Neuronal systems - Bifurcation curves. 


1 Introduction 

The FitzHugh-Nagumo (FHN) oscillator model [Tl|2] is a 
mathematical description of the qualitative features of the 
nervous impulse conduction in the neural cell, namely the 
membrane potential dynamics, the resting bias of the ex¬ 
citable membrane and the existence of an action poten¬ 
tial threshold. The two-dimensional FHN oscillator is de¬ 
scribed by the variables u and re, that represent the poten¬ 
tial of the cell membrane performing the excitability of the 
system, and the recovery variable after an action poten¬ 
tial, respectively. Following Campbell and Waite [3], such 
model is represented by the following differential equations 


du 

dt 

dw 

dt 


= c{w -^u- 

o 

= —(u — abw), 
c 


( 1 ) 


where J(t) is a time-dependent function that represents 
the external stimulus, and a, 6, and c are parameters [U 
|2]. Such system can produce chaotic motion only if an 
external stimulus J(t) is applied [HIS]. Oscillatory chaotic 
behaviors can also be induced considering the coupling of 
FHN oscillators in a network 
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Here, we consider two coupled FHN oscillators with 
two types of coupling, namely unidirectional and bidi¬ 
rectional. Our main interest is to obtain the bifurcation 
curves and consequently in describing the dynamics of the 
variables regarding the coupling strength and other sys¬ 
tem parameters in the parameter-planes. The parameter- 
planes are diagrams when two parameters of the system 
are varied and the other are kept constant. The colors 
codify some quantity that can be computed on the model. 
Usually, this quantity is the Lyapunov exponent [8 | l ^fTQ} 
HB, periods or other invariant quantities iniiii]. 

That procedure allows us to identify regions of periodic 
and chaotic behavior, and recently it is applied in several 
models [HUTMTHUTt] . a general feature is observed in these 
parameter-planes, the existence of shrimp-shaped periodic 
structures embedded in chaotic domains niiiTj. 

In this paper we investigate the dynamics of two cou¬ 
pled FHN oscillators as a function of the type of cou¬ 
pling, without external stimulus. The dynamics is inves¬ 
tigated in the parameter-planes, namely via Lyapunov 
and isoperiodic diagrams and we com¬ 

pare their structural changes as the type of coupling. It 
is shown the limitation of the Lyapunov diagrams to dis¬ 
tinguish the quasi-periodic from periodic oscillations. For 
two FHN oscillators unidirectionally coupled we report pe¬ 
riodicity domains similar to the Arnold tongues. The role 
of the coupling in FHN-networks was studied in recent 
papers but with different approaches of those 





2 


Anderson Hoff et aL: Numerical bifurcation analysis of two coupled FitzHugh-Nagumo oscillators 


studied here. In those papers, the interesting is in the syn¬ 
chronization phenomena, and in our work we carry out a 
systematic dynamical study of the two FHN oscillators 
with respect to the coupling type. 

This paper is organized as follows. In Section the 
network models are introduced, and in the Section^ we 
present the numerical results concerning with each cou¬ 
pling type. In Section we present the conclusions of our 
work. 


2 FitzHugh-Nagumo network models 


In this section we present the FHN-network models for two 
types of coupling between two oscillators, namely unidi¬ 
rectional and bidirectional coupling. 

Following Campbell and Waite [3], for the unidirec¬ 
tional coupling case between two oscillators, the FHN- 
network model is given by 


dxi 

dt 

dyi 

dt 

dX2 

dt 

dm 

dt 


= c{yi +xi- ixf) + 7(a;i - X2), 

= -t{xi -a + byi), 
c 

= c{y2+X2 - ^cl), 

= -^{X2 - a + by2), 
c 


( 2 ) 


where Xi and yi, i = 1 , 2 , represent the voltage across 
the cell membrane, and the recovery state of the resting 
membrane of a neuron, respectively. On the other hand 
a, 6 , and c are parameters, and 7 is the coupling strength 
between the network elements. 

For bidirectional coupling case between two oscillators, 
the FHN-network model can be written as 


dxi , 1 o, , , 

= c{yi -h - -Xi) + 7 (xi - X 2 ), 

-TT = —(^1 -a + byi), 
dt c 

dX 2 / ^ / \ 

= c{y2 +X2- -X2) + 7(^2 - Xi), 

^ = --{X2 -a + by2). 
dt c 


Variables and parameters have the same meaning as in 
system (§. 

In the next section we present and discuss the results 
obtained from the numerical solutions of those equations 
whose behaviors depend on the four above-mentioned pa¬ 
rameters. The results are essentially presented in two- 
dimensional diagrams, using the Lyapunov exponents and 
the periods as the invariant quantities, here namely Lya¬ 
punov and isoperiodic diagrams, respectively, as reported 
previously nans]. Bifurcation curves obtained by numer¬ 
ical continuation method [23] are also presented for both 
network models. 

The Lyapunov diagram plots are obtained for the fol¬ 
lowing parameter combinations: 7 x a with b and c fixed. 


7 x 6 with a and c fixed, and 7 x 0 with a and b fixed, for the 
systems § and (§ . They were constructed using the Lya¬ 
punov exponents numerically calculated for these systems. 
To evaluate the Lyapunov spectra, the Eqs. and @ 
are numerically solved by the Runge-Kutta method with 
fixed time step in 1.0 x 10 “^ and 5.0 x 10 ^ integration 
steps, via the algorithm proposed in Ref. [24] for each pa¬ 
rameter pair discretized in a grid of 500 x 500 values. We 
performed other tests using smaller time steps and the 
results are practically unchanged. Therefore, we obtain 
2.5 X 10^ Lyapunov spectra for each two-dimensional dia¬ 
grams, where each spectrum has the number of values of 
Lyapunov exponents equal to the number of dimensions 
of the network. 

The isoperiodic diagram plots are also computed for 
the two coupled FHN oscillator models. To evaluate the 
period for each time-series with a fixed set of parame¬ 
ters, we use the Runge-Kutta method again, but now with 
variable time step and removed a transient time equal to 
5 X 10^ and 1 x 10^ integration steps to find the period 
obtained by the maxima of the time-series with a preci¬ 
sion of 1 X 10“^. The parameter pairs were discretized in 
a grid of 10 ^ x 10 ^ values. 

To obtain the bifurcation curves in the parameter- 
planes, we used the numerical continuation method. For 
this purpose, the MATCONT package [23] was used, and 
Hopf, saddle-node, Neimark-Sacker, and period-doubling 
curves were obtained unveiling the bifurcation structures 
of the systems Q and The bifurcation curves are 
used extensively to unveil the bifurcation structures of 
dynamical systems |25[l26[l27[l28p29l[3Q] . In one of the first 
papers on this matter m, Barrio et al performed a de¬ 
tailed study of bifurcations in the three-parametric three- 
dimensional Rossler model. Among other results, the au¬ 
thors overlapped bifurcation curves, obtained by numeri¬ 
cal continuation method [23] , on two-dimensional parame¬ 
ter-spaces of the model, and some periodic structures might 
be described by those curves, as also be done recently in 
Ref. [30], for a four-dimensional Chua model. 


3 Numerical results 

We begin our numerical investigation analyzing the para¬ 
meter-planes of systems and (§ in a large range of 
parameter values. The main goal here is to know how is 
the dynamical behavior of the systems regarding the bi¬ 
furcation curves and Lyapunov exponents in those param¬ 
eter values. Based in the results presented in Ref. [3], we 
show in Fig. H the parameter-planes of the systems © 
((a)-(c)) and Q ((d)-(f)) for the largest Lyapunov expo¬ 
nent with bifurcation curves overlapped. The set of initial 
conditions used is (:ri, ^ 1 , ^ 2 , ^ 2 ) = (—0.1, 0.1, 0.5,—0.3), 
and the fixed parameters used are (a, 6 , c) = (0.7, 0.4, 2.0), 
as the case. The cyan, red and magenta curves corre¬ 
spond to Hopf bifurcations and the green curve corre¬ 
sponds to limit-point bifurcation. The white, black and 
yellow regions correspond to equilibrium points, periodic 
and chaotic motions, respectively, and codify the sign of 
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Figure 1. Bifurcation curves in the parameter-planes for two FHN oscillators with unidirectional coupling, system ([2| [top line 
(a)-(c)], and for two FHN oscillators with bidirectional coupling, system § [bottom line (d)-(f)]. (a) 7 x a plane, with b = 0.4 
and c = 2.0. (b) 7 x 6 plane with a = 0.7 and c = 2.0. (c) 7 x c plane with a = 0.7 and b = 0.4. (d) 7 x a plane, with b = 0.4 and 
c = 2.0. (e) 7 X 5 plane with a = 0.7 and c = 2.0. (f) 7 x c plane with a = 0.7 and b = 0.4. The colors white, black, and yellow 
represent the largest Lyapunov exponent and correspond to equilibrium points, periodic and chaotic motion, respectively. The 
arrows indicate small regions that present chaotic motion and are enlarged in Fig. 


the largest Lyapunov exponent of the spectrum. For exam¬ 
ple, white corresponds to negative, black to null and yellow 
to positive largest exponent. In Fig. [^the Hopf bifurca¬ 
tions occur when an equilibrium point (white regions) lose 
its stability and a limit cycle (black regions) born. This 
behavior is corroborated by the largest Lyapunov expo¬ 
nents, white and black regions in the diagrams of Fig. 

As can also be observed, the chaotic behaviors (yellow re¬ 
gions) are small regions in the expanded diagrams and the 
arrows indicate where those regions are localized. 


The equations used to describe the dynamical behav¬ 
ior of the FHN-network with two unidirectional coupled 
oscillators are given by the system Q. For the same set of 
initial conditions used in Fig. we show in Figs, j^a)-(c) 
the Lyapunov diagrams for the following planes: (^ 7 x a 
with b and c fixed, (b) 7 x 6 with a and c fixed, and (c) 7 x c 
with a and b fixed. Such diagrams are amplifications of the 
small regions indicated by the arrows in Figs, [^a)-(c), 
and are constructed with the largest Lyapunov exponent 
of the spectrum codified in colors. For example, in Fig.[^ 
yellow-red-blue colors identify chaotic behavior, since for 
points taken from these regions the largest exponent of 
the spectrum is positive. In Figs, [^d)-(f) we present the 
Lyapunov diagrams for the largest exponents of the spec¬ 
tra for two bidirectionally coupled models, described by 
system Such diagrams are amplifications of the small 
regions indicated by the arrows in Figs. Bd) -(f). As the 
unidirectional case, the Lyapunov spectrum has four ex¬ 
ponents, and the diagrams shown in Figs. |^d)-(f) were 
constructed with the largest one. The construction and 
interpretation of colors follow the same of Figs, [^a)-(c). 


All diagrams in Fig. [^a)-(f) show large periodic re¬ 
gions (regions in black) and chaotic windows (yellow-red- 
blue regions) with periodic structures (in black) embedded 
on them. As far as our knowledge, this feature is commonly 
found in almost every nonlinear dissipative dynamical sys- 

the existence of paths crossing the chaotic regions in which 
the variables and yi have periodic behaviors. Moreover, 
there exist regions in which these variables have chaotic 
oscillations, depending of the parameter combinations. It 
is important to observe that for the single FHN oscillator 
without external forcing, namely a single neuron with¬ 
out external stimulus, system Q, the variables can not 
present chaotic motion. 


The limitations of the Lyapunov diagram shown in 
Fig. I^are related to the distinction between periodic (in 
black) and chaotic (yellow-red-blue) regions, and to ob¬ 
serve the presence of periodic structures, as cuspidal, non- 
cuspidal and shrimp-like structures m. In such dia¬ 
grams is not possible to infer the periods of the struc¬ 
tures and their organization rules, as period-adding cas¬ 
cades and Stern-Brocot trees fT2llT3l[T6| . Fol¬ 

lowing recent works about isoperiodic diagrams [uni 
US], that are parameter-planes using the numbers of spikes 
of the time series (periods), we also constructed these di¬ 
agrams for regions inside the boxes of Figs, [^a) and (c). 

Figure shows the Lyapunov and the isoperiodic di¬ 
agrams for the boxes A and B in Fig. I^a). In these dia¬ 
grams it is clear the limitations of the Lyapunov diagram 
in both cases. The isoperiodic diagram in Fig.j^c) unveils 
hidden regions and pieces of periodic structures that are 
not perceived in Fig. [^a). For example, on the left top in 
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Figure 2. Lyapunov diagrams for the largest exponent of the spectra codified by colors, as the color bar at the right side, for 
two FHN oscillators with unidirectional coupling, system § [top line (a)-(c)], and for two FHN oscillators with bidirectional 
coupling, system ^ [bottom line (d)-(f)]. The diagrams are amplifications of small regions located by the arrows in Fig. 
Periodic/quasi-periodic domains are represented by black colors, and chaotic ones are yellow-red-blue colors. The boxes in panels 
(a) and (c) delimit the regions magnified in Figs. and respectively. 


Fig-i a) we observe a large black region with branches 
of periodic structures coming out from it and extending 
over the chaotic region (yellow-red-blue region), however, 
in Fig. I^c) we observe in this same region a combination 
of quasi-periodic behavior (black color) with branches of 
periodic structures embedded on it. Moreover, the orga¬ 
nization pattern of these periodic structures delimits the 
quasi-periodic region, upper black region in Fig.l^c), with 
the chaotic region, lower black region in Fig. [^CL This or¬ 
ganization pattern is similar to the Arnold’s tongues [32l 
l33] . and its periodicity is similar to branches of Stern- 
Brocot tree observed in a variety of continuous-time sys¬ 
tems [T2j[T3l[T6l|34] . Basically, in Fig. [^c) the organiza¬ 
tion rule of the periods can be summarized as follows. 
The primary structures follow a period-adding cascade: 
7^6^5^4^3^2. The sum of the periods of two 
consecutive primary structures is equal of the period of 
secondary structures that are between these two primary 
structures. For example: 7 + 6 ^ 13, 6 + 5^ 11,5 + 4^ 
9,4 + 3^ 7, 3 + 2^5 (see Fig. |^c)). The sum of the 
period of a primary structure with the consecutive sec¬ 
ondary structure is equal of the period of tertiary struc¬ 
ture between these primary and secondary structures. For 
example: 7 T 13 —y 20, 6 T 13 —y 19, 6 T 11 —y 17, 5 T 11 —^ 
16,5 T 9 — y 14,4 T 9 —y 13,4 T 7 — y 11,3 T 7 — y 10,3 T 5 — y 
8 ,2 T 5 —y 7. 

To corroborate the differences between the three dis¬ 
tinct behaviors presented in Figs. |^a) and (c), we show 
in Fig. three attractors in three different points of the 
Lyapunov diagram of Fig. [^a) (see the points I, II, and 
III). In each plot we also show the respective Lyapunov 
spectrum (see the insets in Fig. [^. 

FiguresJ^b) and (d) show the amplification of the box 
B in Fig. I^a), and another pattern formation of peri¬ 
odic structures can be visualized. Here, the structures self¬ 
organizing in a circular layers pattern, in which each layer 


is composed by two connected structure with same period. 
A period-adding cascade occurs in direction of the inner 
layers. 

Figures l^b) and (c) show the two other planes, (b) 
7 x 6 , and (^7 x c for the largest Lyapunov exponents of 
the system Q. In Fig.j^b) the similarities with the 7 x a 
plane (see Fig.j^a)), are evident with the two types of or¬ 
ganization patterns of structures shown in Fig.[^ In 7 x c 
plane (see Fig. [^c)), we observe the same periodic struc¬ 
tures similar to the Arnold’s tongues observed previously 
in the 7 x a plane (see Fig.|^. For example, Figs.j^a) and 
(b) are amplifications of the Lyapunov and isoperiodic di¬ 
agrams, respectively, for the box in Fig. [^c). Bifurcation 
curves are also overlapped in the Lyapunov diagram of 
Fig. E a), where the cyan and blue curves are saddle-node 
bifurcations, and the green curves are Neimark-Sacker bi¬ 
furcation curves. It is clear to observe the same organiza¬ 
tion pattern presented in Fig. |^a) and (c), quasi-periodic 
regions (lower left and right portions of the diagram) with 
branches of periodic structures embedded on it, and other 
branches of these structures embedded on chaotic regions 
(higher portion of the diagram). The same organization 
rule of periods of these structures was also observed. To 
corroborate this statement, in Fig. we show amplifica¬ 
tions of Fig. I^and we observe similarities with Figs, [^a) 
and (c). In Figs. |^b) and (d) the same organization rule 
of the Arnold tongues in Figs.j^a) and (c), or in Figs.J^a) 
and (c), is observed for the periodic structures of Figs, [^b) 
and (d). 

Another remarkable result is the agreement between 
the bifurcation curves in Fig.[^a) and the periodic struc¬ 
tures. In addition, there exist a connection between the 
lower left Arnold tongues and the central periodic struc¬ 
tures. Initially, in Fig. I^a) it is clear to see two saddle- 
node bifurcation curves (cyan curves) bordering two cen¬ 
tral periodic structures (the central yellow and red peri- 
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Figure 3. Unidirectional coupling case: (a) and (b) Lyapunov diagrams for the largest exponents, and (c) and (d) isoperiodic 
diagrams, both for the boxes A and B, respectively, in Fig. |^a). In (a) and (b) the black-yellow-red colors codify the largest 
Lyapunov exponent, and in (c) and (d) the top color bars codify the periods. Black corresponds to quasi-periodic/chaotic 
behaviors or periods greater than 20. 



Figure 4. Attractors of three distinct points in Fig.|^a), showing (a) Torus-2 (point I), (b) periodic (point II), and (c) chaotic 
(point III) behaviors. The Lyapunov spectrum for each attractor is also shown. 


odic structures in Fig.j^b)) at their lower portions. Those 
saddle-node bifurcations occur when, as the control pa¬ 
rameters are varied, the chaotic attractors become unsta¬ 
ble and born stable periodic attractors. Those features are 
corroborated by the isoperiodic diagram in Fig. [^b). Fol¬ 
lowing the second cyan curve (the curve that border the 
red structure in Fig. [^b)) from right to left, we see that 
after to border the lower portion of the periodic structure 
this curve reach the respective left Arnold tongue of same 
period (the red Arnold tongue in Fig.l^b)). Due to multi¬ 
stability the yellow Arnold tongue was hidden by the lower 


blue periodic region (see Fig. |^b)). Moreover, due to nu¬ 
merical precision to obtain those bifurcation curves, how 
smaller is the periodic structure, more difficult is to ob¬ 
tain the curves. Therefore, for period-4 structure and on, 
we have not accuracy to numerically obtain their saddle- 
node bifurcation curves. However, our results show that 
the formation of the central set of periodic structures is 
connected with the formation of the left set of Arnold 
tongues. 

The blue saddle-node bifurcation curve in Fig. [^a) 
delimits the blue periodic region in the bottom of the 
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Figure 5. Unidirectional coupling case: (a) Lyapunov diagram for the largest exponents with bifurcation curves overlapped, 
and (b) isoperiodic diagram, both enlarged of Fig. |^c). In (a) the cyan and blue curves are saddle-node bifurcation curves, 
and the green curves are Neimark-Sacker bifurcation curves. In (b) the top color bar codifies the periods. Black corresponds to 
quasi-periodic/chaotic behaviors or periods greater than 20. 


isoperiodic diagram of Fig. ^b). The horizontal green 
Neimark-Sacker curve in Fig.^a) (in the bottom) shows 
the line where born the Arnold tongues, corroborated by 
the isoperiodic diagram in Fig. I^b), where in the bot¬ 
tom of the diagram we see the branches of the Arnold 
tongues in the quasiperiodic domains (lower black regions 
in Fig.l^b)). The diagonal green Neimark-Sacker curve in 
Fig. [5[a) has a strict connection with the multistability 
presented by the system (§ and with a family of Torus 
that emerges very close to the curve. To show this feature, 
in Figs.J^a) and (b) we present in more details a zoom 
of Fig. Imb), with the green Neimark-Sacker and cyan 
saddle-node bifurcation curves for two different initial con¬ 
ditions: (a) (xi, ^ 1 , X 2 , ^ 2 ) = (—0.1, 0.1, 0.5, —0.3), and (b) 
(xi, ^ 1 , X 2 , ^ 2 ) = (—0.1, 0.5, 0.1,—0.3). By contrasting the 
diagrams in (a) and (b) we see the effect of the initial con¬ 
dition and the role of the green Neimark-Sacker curve, for 
different initial conditions the black periodic border in¬ 
creases until near the limit of the green curve destroying 
chaotic regions and periodic structures, and very close to 
the green curve a Tori domain emerge. The Tori domain 
is observed in Fig.j^c) by plotting the Lyapunov diagram 
with the second largest exponent. In this case, the periodic 
domains are in white. Tori (quasiperiodicity) and chaotic 
ones are in black. The Tori domain is the black strip just 
above the Neimark-Sacker bifurcation curve (green curve), 
once that this strip is black in the diagrams of Figs. lilb) 
and (c). To corroborate, in Fig.j^d) we present the isope¬ 
riodic diagram of Fig. [^b). In that diagram, near to the 
border of the blue region emerge Arnold tongues and a 
Tori domain close to the Neimark-Sacker curve, delimit¬ 
ing the boundary between the Tori and chaotic domains. 

In Figs. I^a) and (b) we show two planes of the four¬ 
dimensional basin of attraction of the point P in Figs.l^a) 
and (b). The basins are constructed using the largest Lya¬ 
punov exponent and associating black and red colors for 
periodic and chaotic attractors, respectively. The initial 
conditions in Fig. |^a) are (^ 2 ,^ 2 ) = (0.5,—0.3), and in 


Fig. [^b) (xi,^i) = (—0.1,0.!). It is easy to observe the 
existence of large periodic domains (black regions) in the 
basins of attraction. 

Regarding the two coupled FHN oscillators with bidi¬ 
rectional coupling. Figs, [^d)-(f), the dynamics does not 
display the richness of details shown in Figs. ia)-{c), 
for the two coupled FHN oscillators with unidirectional 
coupling. Comparing the dynamics, presented in the Lya¬ 
punov diagrams of Fig. for both coupling, it is possi¬ 
ble to conclude that the unidirectional coupling presents 
a more subtlety in the dynamics concerning the presence 
of self-organized periodic structures embedded in quasi¬ 
periodic and chaotic regions in the three planes of param¬ 
eters, namely 7 x a, 7 x 6 , and 7 x 0 . 

4 Conclusions 

In this paper we have investigated the dynamics of FHN 
oscillators composed by two coupled models, using uni¬ 
directional and bidirectional coupling. Each oscillator is 
modeled by a set of two autonomous nonlinear first-order 
ordinary differential equations that describes the dynam¬ 
ics of a nerve impulse through the neuronal membrane, 
namely FitzHugh-Nagumo model. Our main goal is the 
investigation of the influence of the coupling strength be¬ 
tween the oscillators with the type of coupling. For this 
purpose, we constructed the Lyapunov and isoperiodic di¬ 
agrams, which are parameter-planes with the Lyapunov 
exponents and periods codified by colors, respectively, for 
the two cases of the FHN-networks. By contrasting those 
diagrams in each case, we observe the changes in the dy¬ 
namical behaviors of these models. We showed that 7 
also has an important role in the dynamics of the FHN- 
networks. In addition, we have obtained the bifurcation 
curves by numerical continuation method, and overlapped 
those curves in the Lyapunov diagrams, unveiling the ex¬ 
istence of a connection between different types of periodic 
structures through those bifurcation curves. 
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Figure 6. Unidirectional coupling case: (a) and (b) Lyapunov diagrams for the largest exponents, and (c) and (d) isoperiodic 
diagrams, both for the box B and C, respectively, of Fig.[^ In (a) and (b) the colors follow the same of Fig.[^ and in (c) and 
(d) the top color bars codify the periods. Black corresponds to quasi-periodic/chaotic behaviors or periods greater than 20. 


Regarding the two coupled FHN oscillators for the uni¬ 
directional coupling, Le., the coupling only occurs in one 
of the four variables {xi^yi]i = 1,2), namely xi in sys¬ 
tem (§. The dynamics for this case in the Lyapunov and 
isoperiodic diagrams presents the coexistence of periodic, 
quasi-periodic, and chaotic oscillations, and the presence 
of stability domains (periodic structures) embedded in 
quasi-periodic and chaotic domains. The dynamical be¬ 
havior presented in this system resembles with the dynam¬ 
ics of systems that show the presence of Arnold tongues. 
The organization rule of periods, shown in Figs.j^b) and 
l^c), seems to be a branch of the Stern-Brocot tree previ¬ 
ously reported in a three-dimensional version of the FHN 
model [12]. A circular organization of periodic structures 
in layers was also observed in a portion of the parameter- 
plane (Figs.j^b) and (d)). 

For the bidirectional case, the coupling occurs in two 
of the four variables {xi^yi]i = 1,2), namely xi and X 2 
in system In this case, the dynamics change drasti¬ 
cally with the disappearance of the quasi-periodic domain 
and of the Arnold tongue-like structures, surviving the 
chaotic domains and some periodic structures (see Fig.[^. 
Among some examples of discrete- and continuous-time 
systems that present Arnold tongues in the parameter- 
planes forced oscillators show a large spectrum 

of applications and Arnold tongues are abundantly ob¬ 


served in those systems. In our work we characterize the 
Arnold tongues in the parameter-planes of the unidirec¬ 
tional coupling case that, indeed, can be interpreted as 
a forced oscillator, being (^ 2 ,^ 2 ) the master and (xi^yi) 
the slave in system (§. The primary Arnold tongues and 
periodic structures shown in the isoperiodic diagrams of 
Figs. [^c), 1^1^ and[^ are organized by period-adding cas¬ 
cades with period-1 being the factor of the adding. An in¬ 
teresting verification is that the master oscillates in cicle-1 
driving the slave with a period-1 oscillations, which is the 
factor of the period-adding cascades of the slave. 

For networks of two models, and for both coupling 
cases, the Lyapunov diagrams are slightly different over 
wide regions (see Fig. but a general feature was ob¬ 
served in some regions of these diagrams, namely the ex¬ 
istence of periodic structures embedded in chaotic regions. 
These sets of periodic structures are presented in a wide 
range of nonlinear systems j9lll5[ll6[ll7j . An exception is in 
high-dimensional systems with more than three-dimensions, 
where hyperchaotic behaviors can occur. In these systems, 
on the hyperchaotic regions there are no periodic struc¬ 
tures or the structures are malformed or shapeless [TS] , 

Relating the study presented here with the synchro¬ 
nization regimes in the two coupled FHN systems, it was 
also observed that, for both systems (unidirectional and 
bidirectional coupling) initializing in the same initial con- 
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Figure 7. Unidirectional coupling case: (a) and (b) Lyapunov diagrams for the 7 x c plane for the largest exponents for two 
distinct initial conditions, and in (c) Lyapunov diagram for the second largest exponent for the same initial condition in (b). 
Cyan and green curves are saddle-node and Neimark-Sacker bifurcations, respectively. 


(a) 

4- 

(b) 
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- 4 - 
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Figure 8 . (a) and (b) planes of the four-dimensional basin of attractions for the point P in Figs.[^a) or (b). Black color is the 
basin of periodic attractors, and red one for chaotic attractors. 


ditions, lead to synchronized states. On the other hand, 
initializing both systems (unidirectional and bidirectional 
coupling) in different initial conditions, lead to non-syn¬ 
chronized states in the periodic domains of the parameter- 
planes of Fig. [2 Moreover, chaotic synchronization was 
not observed in the chaotic domains covered in our study. 
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